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Abstract. Recent work on analytical solutions to the MSW equations for three neutrino 
flavours is reviewed, with emphasis on the exponential density. Application to a particular 
mass matrix, proposed by Lehmann, Newton and Wu, is also discussed. Within this model, 
the experimental data allow a determination of the three neutrino masses. They are found to 
be 0.002-0.004, 0.01 and 0.05 eV. 



1. Introduction 



We here review some recent results on analytical solutions of the Mikheyev-Smimov- 
Wolfenstein (MSW) effect [[I]] for the propagation of three neutrino flavours. Such analytic 
results have been obtained for both the exponential density [||, ^ and the linear density [0]. 

In the case of an exponential electron density, the three neutrino wave functions can be 
expressed in terms of generalized hypergeometric functions, 0- linear density, 

the solutions can be expressed in terms of a Fourier transform In the case of two flavours, 
these reduce to parabolic cylinder functions or confluent hypergeometric functions, but of 
different parameters and arguments [^, ^ . 

A particular neutrino mass matrix, originally proposed by Lehmann, Newton and Wu 
(LNW) for quarks [0] is also reviewed [|], ^ [1^]. Within this model, the current data on 
atmospheric and solar neutrinos permit a determination of the neutrino masses. 



2. Exponential density 



In a medium where the electron neutrino (described by <^i{r)) interacts differently from the 
others, the propagation is given by the equation 
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The mass matrix is real and symmetric, Mjj 



M^- 

'J 



{M~)-j, but otherwise arbitrary. (The 



LNW mass matrix will be discussed later.) Furthermore, D{r) = \/2GY;Ne{r), with Gp the 
Fermi constant and Ne{r) the solar electron density. 

For the sun, the density [ [TT| ] is well approximated by an exponential, 

Ne{r)=Ne{0)e-'''^\ Tq^OAxEq. (2) 

Introducing a new radial variable: w = r/rp + Mq, and performing a rotation on the second 
and third components, Eq. (|l]) can be written as 
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The eigenvalues of the 3x3 matrix 
«i X2 

X3 «3 

are denoted /ij, and /ig, they are the squares of the neutrino masses multiplied by r^/ {2p). 
Together with cOj and (O2 they control the evolution of the i//;. 

Introducing now the variable z = ie^", the solutions to Eq. (Q) can be expressed in terms 
of solutions to the third-order ordinary differential equation 
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namely generalized hypergeometric functions []3|, p^]: 
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The solutions to Eq. (|^) are thus 

v/, = qvA(i)+qv/(2)^C3VA(3), 

where the constants C- are determined by the boundary conditions: i/Zj = 1, i/Zj = V^3 = at 
the center of the sun, r = 0. 

The 2^2 functions can be defined in terms of the series expansions 
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where (a)^ is a Pochhammer symbol, but (somewhat complicated) asymptotic approxima- 
tions are actually more useful for numerical work |^]. 



3. The LNW mass matrix 

For quarks, it was found [0] that a particular, simple texture for the d (d, s, b) and u (u, c, t) 
quark mass matrices leads to a good description of the CKM matrix [[T3[]. The mass matrix is 
assumed to have the form 



M 



d 
deb 
b a 



(9) 



with b^ = 8c^. This is known as a two-zero texture, but differs from those normally considered 
(see, e.g., |T5]]) in the additional relation between b and c. The eigenvalues are given by 
mj, and m^, with mj < m^. 
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Actually, in the case of quarks, since there is CP violation, a complex CKM matrix 
is obtained by replacing the parameter d in @ by ±?J for u quarks (such that M remains 
Hermitian). The CKM matrix becomes 

^CKM = ^''(")diag(-/, 1, l)R{d), (10) 

The Jarlskog determinant [ [T^ ] thus obtained is 7 = 2.6 x 10^, in good agreement with data. 

This same mass matrix has been applied to the case of three neutrinos ^, [T^], and 
rather good fits to the atmospheric [jFTp and solar [jTS], |T9|, ^ ^ ^ neutrino data 

have been obtained. 

The matrix M is diagonalized, whereby M = RM^^^^R^. For this purpose, the following 
notation is convenient 

5j =m3 — m2 + mj, 
= a + c 

— ^2 = m^m2 — m^m^ + m2m^, 
ac 



Sc^+d^ 



ad 



(11) 



Eliminating d and c, the resulting cubic equation for the parameter a can be written as 

9a^-llS^a^ + {SSj + S2)a-S^ = 0. (12) 

A physical solution requires a real and positive. This is equivalent to having three real 
solutions for a. One of these is negative and two are positive. At any point inside the allowed 
domain in the {m^/m-^)-{m2/m-^) plane (See Fig. 1), there are thus two allowed solutions, 
denoted Solutions 1 and 2. 

This diagonalization has to be carried out for both neutrinos and charged leptons, in order 
to obtain the neutrino mixing matrix [ [T0| ] 
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(13) 



where (e = 1 or z, depending on whether or not there is CP violation) 

t/ = (^CKM)^=^^Wdiag(e,l,l)7?(v) 
relates the neutrino mass eigenstates to the flavor states of charged-current interactions: 

I V.) = ^,1 1 Vi ) + U^2 1 ^2) + ^.3 1 ^3) ' etc. 

There are four possible solutions to Eq. ([T2|), two for the neutrino sector, and two for 
the charged lepton sector. Furthermore, the model has two signs (denoted "parities") to 
be specified: b and d can each be either positive or negative (unless d is imaginary [|T0|]). 
However, only the product of the parities" in the neutrino and charged lepton sector matters, 
and similarly for the "J parities", so we may put both "parities" of the charged lepton sector 
to +1. 



(14) 



(15) 
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4. Fits to data 



Let us consider first the atmospheric neutrino data. The Super-Kamiokande results [ [r7| ] give 
Am^ ~ (2 — 3) X 10^-^ eV, with sin^(20) ~ 1. The survival of muon neutrinos is given by 
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(16) 



where U is the neutrino mixing matrix of Eqs. ( [T3[ ) and ([T4]). In the limit of Am^^t/Ap <^ 1 
this simplifies, and invoking further unitarity, one finds 

P,,_,,.fn ~ AUf.^Ut. sin" ( — ^ 1 , (17) 



which suggests that one needs \U^^U^^\ = ^(1). This can be achieved within the model (for 
both Solutions 1 and 2), for <^m^, with also small compared with m^. Furthermore, the 
data suggest that the scale mg must be such that m^ — (2 — 3) x \0^^ eV, or = ^(0.05 eV). 

Fits to atmospheric data confirm this qualitative analysis. Invoking also the solar CI 
[^, Ga 0, Super-Kamiokande and SNO neutrino data, one finds that 

Solution 2 for the neutrino sector. Solution 1 for the charged lepton sector and "Z? parity"=— 1 
give good fits for <C m^, with /nj also small as compared with 7^3. A determined from 
these different atmospheric and solar survival probabilities leads to good fits (see Fig. 1) with 
about 0.052 eV, about 0.01 eV, and ~ 0.002-0.004 eV g The "J parity" is 
unimportant. 



m, [eV] 1 




m, [eV] '= 




ution (2,1 ) — 
m3 = 0,052 eV 



m, [eV] 1 



0,04 
m, [eVl 




Figure 1. Fits to the atmospheric and solar neutrino data for Solution 2 (neutrino sector) and 1 (charged 
leptons) and both "paiities" negative. Contours are given at = 5, 10, 15, 20, 25. Sum (right panel) 
also at 30, 35, 40, 45, 50. The best-fit point (right panel, at the lower left), is marked "low". 

The effect of the CHOOZ data is to disfavor the region m2 = ^{m^), where \U^^ \ is 
large. But this region is already disfavored by the solar neutrino data, as is seen in the middle 
panel of Fig. 1. In terms of the more conventional two-flavour analyses for the solar-neutrino 
sector, these fits roughly correspond to the large-mixing-angle solution. 
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5. Summary 

Analytic results for the solutions to the MSW equations for three neutrino flavours are very 
valuable for a fast scanning over the parameters of some given model for the mass matrix. 

The LNW mass matrix is a very constrained model that in the quark sector describes the 
CKM matrix, and in the neutrino sector gives the mixing in terms of the mass eigenvalues. 
Applied to the neutrino data, the model gives a very good fit. 

The solar neutrino data has also been studied within the same model, using numerical 
integration methods (no 2^2'^) Wh- An additional fit was then found at ~ 2.8 x 10^^ eV, 
corresponding to the small-mixing-angle solutions. However, this point is disfavoured by the 
atmospheric neutrino data, and by the flat electron recoil spectrum [|l9p. 
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